Topologically Protected Quantum State Transfer in a Chiral Spin Liquid 
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Topology plays a central role in ensuring the robustness of a wide variety of physical phenomena. 
Notable examples range from the robust current carrying edge states associated with the quantum 
Hall and the quantum spin Hall effects to proposals involving topologically protected quantum 
memory and quantum logic operations. Here, we propose and analyze a topologically protected 
channel for the transfer of quantum states between remote quantum nodes. In our approach, state 
transfer is mediated by the edge mode of a chiral spin liquid. We demonstrate that the proposed 
method is intrinsically robust to realistic imperfections associated with disorder and decoherence. 
Possible experimental implementations and applications to the detection and characterization of 
spin liquid phases are discussed. 



The decoherence of both quantum states and quantum 
channels represents a major hurdle in the quest for the 
realization of scalable quantum devices [1, 2]. Several 
avenues are currently being explored to address these 
important challenges. For example, quantum repeater 
protocols are expected to improve the fidelity of quan- 
tum state transfer - the fundamental building block of 
quantum communication [3, 4]. Similarly, quantum error 
correction can significantly extend the lifetime of quan- 
tum memories and suppress the errors associated with 
quantum logic operations [5, 6]. The practical realiza- 
tion of these technologies, however, requires a high level 
of quantum control that is as yet, not experimentally ac- 
cessible. An alternative paradigm to achieving protected 
quantum states is provided by topology; indeed, if such 
states can be stored in the topological degrees of freedom 
of certain exotic states of matter, they become intrinsi- 
cally robust against local noise [7-12]. 

The implementation of robust long-lived quantum 
memories can also be achieved by encoding quantum bits 
in appropriately chosen physical degrees of freedom. For 
example, the natural isolation of nuclear spins immunizes 
them from the environment and makes them an excep- 
tional candidate for the storage of quantum information 
[13-18]. Such solid-state spin qubits can be locally cou- 
pled with high fidelities, enabling the realization of few- 
bit quantum registers [14, 16]. However, spatially remote 
registers interact extremely weakly; thus, in this context, 
the challenge of scalability is shifted to the development 
of quantum channels capable of connecting remote regis- 
ters in a robust and noise-free fashion [19-21]. 

This article describes a novel approach to the realiza- 
tion of intrinsically robust quantum channels and exploits 
topological protection to enable high-fidelity quantum in- 
formation transport. We envision quantum state transfer 
between remote spin registers to be mediated by a 2D sys- 
tem composed of interacting spins. Specifically, the spin 



system is tuned into a gapped chiral spin liquid phase, 
which harbors a fermionic edge mode. The prototype 
of this specific chiral spin liquid is the gapped B phase 
(CSLB) of the Kitaev honeycomb model [22]. Although 
such a phase is best known for its non-Abelian vortex 
excitations, here, by operating at finite temperatures 
below the gap, we make use of its Majorana fermionic 
edge mode as a topologically protected quantum channel. 
Moreover, we discuss possible applications of our proto- 
col for the spectroscopic characterization of spin liquid 
states [23]. 



APPROACH TO TOPOLOGICALLY 
PROTECTED STATE TRANSFER 

Our approach to quantum state transfer is schemati- 
cally illustrated in Fig. 1. Quantum information is en- 
coded in a two-qubit spin register, with each qubit ca- 
pable of being individually manipulated. The register is 
coupled to the edge of a two-dimensional spin droplet, 
whose elements we assume cannot be individually ad- 
dressed but can be globally "engineered" to create a spin 
liquid state in the CSLB phase. The transfer protocol 
proceeds by mapping the quantum information stored in 
the left-hand spin-register onto the chiral edge mode of 
the droplet. The resulting wavepacket traverses the edge 
before retrieval at the remote register. 

A distinct feature of our protocol, as compared with 
previous approaches [20, 21, 24-26], is the fundamen- 
tal robustness of the quantum channel. The chiral na- 
ture of the fermionic edge mode ensures that destruc- 
tive backscattering during state transfer is highly sup- 
pressed; moreover, the characteristic (linear) dispersion 
of the edge-mode ensures that wave packet distortion is 
minimized. Finally, we demonstrate that our approach 
is remarkably insensitive to disorder and decoherence af- 
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FIG. 1. Schematic representation of topologically pro- 
tected state transfer | The grey droplet represents a 2D 
array of interacting spins tuned into the CSLB phase. Quan- 
tum spin-registers composed of a transfer qubit (green) and a 
memory qubit (gold) are arranged around the edge of the 2D 
droplet and coupling between them occurs through the chiral 
edge mode. (1) By mapping the quantum information onto 
a fermionic wave-packet (blue) traveling along the edge, the 
quantum state can be transferred to a remote register. The 
wavepacket travels only in the direction of the blue arrow; this 
chirality prevents mode localization and destructive backscat- 
tering. At a specified time at the remote register location, 
the coupling is turned on and the wavepacket is captured (2). 
Given an ancillary memory qubit and local register manip- 
ulations, a two-qubit gate (3) can be performed before the 
quantum state is transferred back to the original register and 
stored (4-5). This allows for universal computation between 
the memory qubits of spatially separated registers. 



fecting both the bulk and edge of the droplet. Although 
any spin system with a stable CSLB-like phase can po- 
tentially mediate topologically protected state transfer 
(TPST), to illustrate the microscopic mechanism respon- 
sible for such state transfer, we turn initially to a particu- 
lar model and will later generalize our analysis to include 
the effects of disorder, additional interactions, and deco- 
herence. 



TPST ON THE DECORATED HONEYCOMB 

We now consider a specific exactly solvable spin- 
1/2 model which supports robust TPST [27]. Within 
this Yao-Kivelson model, the spins are situated on a 
triangular-decorated honeycomb lattice as depicted in 
Fig. 2 [27]. The associated Hamiltonian naturally gen- 
eralizes the Kitaev model [9] and features a chiral spin 
liquid ground state (CSLB phase), 

x,x' y ,y' z,z' 

links links links 



where a are Pauli spin operators {h = 1). The model 
may be solved by introducing four Major ana opera- 
tors, {7^,7^,7^,7^} for each spin, as shown schemat- 
ically in Fig. 2a and by representing the spin algebra 
as: cr^ = ^7^7^ cr^ = 27^7^ = ij^j^ [22, 27]. The 
Majorana operators are Hermitian and satisfy the stan- 
dard anticommutation relation {7^7^} = ^Sim- The 
Hilbert space associated with the physical spin is a two- 
dimensional subspace of the extended four-dimensional 
Majorana Hilbert space; thus, we must impose the gauge 
projection, P = where D = j^j'^j^j^ [22]. 

Transforming to Majorana operators results in the ex- 
tended Hamiltonian 

= ^T^^ijihl (2) 

where Uij = ^7f 7^ {a depends on the type of ij-link) 

for ij connected and zero otherwise; these Uij correspond 
to the boxed Majorana pairs illustrated in Fig. 2a. Re- 
markably each Uij commutes with the Hamiltonian and 
with all other Ui^rm implying that the extended Hilbert 
space can be divided into sectors corresponding to static 
choices of {Uij = ±1} [22, 27]. 

The choice of {Uij} yields a Hamiltonian which is 
quadratic in the 7^ Majorana operators; from the per- 
spective of these Majoranas, Uij is a static background 
Z2 gauge field. The physical states are sensitive only to 
the flux of the gauge field, w{p) = Hijeap where p 
represents a plaquette, dp is its boundary and ij is ori- 
ented according to the arrows in Fig. 2b [27]. For any 
link with Uij = +1, this orientation can also be inter- 
preted as the direction in which a 7^ Majorana hops in 
order to accumulate a 7r/2 phase. The ground state flux 
sector of the model has w{p) = +1 for all plaquettes, 
corresponding to tt phase around the dodecagonal pla- 
quettes and 7r/2 phase around the triangular plaquettes, 
as shown in Fig. 2b. The 7r/2 phase around the trian- 
gular plaquettes indicates the breaking of time-reversal 
symmetry necessary for a chiral ground state. Alternate 
flux sectors contain plaquettes with vortex excitations 
defined hy w{p) = —1. In general, such vortices are en- 
ergetically gapped by A^, but the energy and dynamics 
of vortices near the edge are controlled by the details of 
the boundary. 

In each flux sector, the associated Majorana Hamilto- 
nian can be diagonalized through a unitary transforma- 
tion Q such that i^^ij Qk,i{'^Uij)Ql^ ^ = Skk'^k, yield- 
ing = ^ ^kclck, where Ck = ^ Qkjlj^ 
N is the number of spins on the lattice, and the index 
k is ordered according to energy. Owing to particle-hole 
symmetry, = —^-k and = c_/c; thus, by restricting 
to /c > 0, 

k>0 

where Ck and satisfy Dirac anticommutation rela- 
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FIG. 2. Coupling between spin-register and the 
droplet edge | (a) Schematic representation of the general- 
ized Kitaev Hamiltonian on the decorated honeycomb lattice. 
Spins are represented by four Majorana operators; spin-spin 
interactions become products of the four Majoranas living on 
each link. Boxed spins correspond to the Uij operators which 
determine the effective hopping associated with the 70 Majo- 
ranas. Coupling (dashed line) between the quantum register 
and the 2D droplet can occur at any edge vertex with an un- 
paired Majorana; the Majorana flavor determines the form 
of the spin-spin interaction that introduces the desired addi- 
tional hopping, (b) Schematic representation of the ground 
state flux configuration. Each arrow represents a Majorana 
hopping of i, yielding tt phase (oriented CCW) around the 
dodecagonal plaquettes and 7r/2 phase around the triangular 
plaquettes. Vortex excitations (circular arrows) correspond 
to the flipping of a Uij link (red arrow and link) , which yields 
uj{p) = — 1 in the two adjacent plaquettes. Quantum registers 
corresponding to a transfer qubit (green) and a memory qubit 
(gold) are shown coupled into the droplet (injection points) 
at two different dangling edge spins. 



tions. Diagonalization of the ground state flux sector on 
a cylinder yields three bulk fermion bands, energetically 
gapped by A5, as shown in Fig. 3 [27]. At the edge, the 
fermionic quasiparticles form gapless chiral modes which 
are guaranteed by the nontrivial Chern number of the 
bulk fermion bands. 



Spin-Register Coupling to a Chiral Edge 

We now consider the addition of spin qubits, which can 
be individually manipulated and read out, to the edge of 
the exactly solved model with open boundary conditions 
[14, 16]. Each edge spin with coordination two has an 
uncoupled Majorana operator, which we term dangling 
as depicted in Fig. 2a. We can extend the definition 
of vortices to include the dangling plaquettes defined by 
the Uij links between dangling Majoranas, as shown by 
the red rectangle in Fig. 2a. These dangling vortices are 
completely decoupled from the fermions and lead to a 
large degeneracy of the model. However, generic per- 
turbations will lift this degeneracy by gapping out these 
dangling vortex states; in this situation, as we later de- 
scribe, the control of dangling vortices at the injection 
point will become important. 

To illustrate TP ST, we consider the full Hamiltonian 
Ht = Hq + Hint where Hint characterizes the coupling 
between the two spin-registers (termed L and R) and 
dangling spins at the edge of the droplet (Fig. 2b), 



H,„t = ~i^L + <^r) + gLCT^a + 9R^>1 (4) 



Here, A 5 is the splitting of the register states (e.g. by an 
applied field), {3^7] are chosen to respect the interaction 
symmetry at the injection points, and Qli Qr represent 
the interaction strength between the registers and the 
injection spins (a and h) as shown in Fig. 2. Transform- 
ing to Majorana operators yields Hint = "^(^TlTl + 
hRlR)^9Llhhha^gR7}i7l7hb^ where, without loss 
of generality, we have chosen a a^a^ register-edge inter- 
action. 

The existence of a dangling Majorana at the droplet 
edge is critical to enable spin-edge coupling. At the injec- 
tion points (Fig. 2b), the register-edge coupling of equa- 
tion (4), not only creates a fermionic excitation, but also 
introduces a dangling vortex (by flipping the Uij cor- 
responding to the adjacent dangling plaquette). Thus, 
in order to exploit the chiral fermion mode to trans- 
port spin-based quantum information, we will need to 
control the injection point (see Supplementary Informa- 
tion for additional details). Imperfections in such con- 
trol will result in the spin-register coupling to additional 
nearby spins. However, since the nearest spins surround- 
ing the injection point will not contain dangling Ma- 
joranas, these additional interactions will naturally gap 
out. 

Even in the presence of the additional interactions pre- 
scribed in equation (4), since UL,a and Ur^ij are con- 
served, the model remains exactly solvable. Expressed 
in terms of the eigenmodes of the unperturbed Hamilto- 
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nian in the ground state flux sector, 



Ht = H^H, 



int = X]e/e(4c/e - ^) 

k>0 



•A5(clcL-^) + A5(4c^-^) 

k k 
k k 



(5) 



where we have deflned ^ 



Z7^'^) and 



CL,i? = 1/2(7^'^ + ^73^'^) 5 ^^^is language, the spin 
state of the L(R) qubit is encoded in the occupation of 
the L(R) fermion mode. The flrst term of the Hamilto- 
nian characterizes the modes of the 2D droplet, the sec- 
ond and third term characterize the splitting associated 
with the spin-registers, while the flnal two terms capture 
the coupling between the registers and the dangling edge 
spins. This Hamiltonian acts in the extended fermionic 
Hilbert space and returning to physical spin states re- 
quires gauge projection (see Supplementary Information 
for details). 



Topologically Protected State Transfer in the Dot and 
Droplet Regime 

The coupling between the register and the chiral edge 
mode can be analyzed in two distinct regimes: 1) the 
mesoscopic dot regime and 2) the macroscopic droplet 
regime. The distinction between these two regimes is 
best understood from a perspective of resolvability; in 
the dot regime, we consider the coupling to a small flnite- 
size system, enabling energy resolution of the individual 
chiral edge modes. Thus, TPST is mediated by a single 
fermionic eigenmode of the system [21]. Meanwhile, in 
the droplet regime, we consider the coupling to a larger 
system, in which energy resolution at the single mode 
level would be extremely difficult. In this regime, we en- 
code the spin register's quantum information in a travel- 
ing fermionic wave-packet. 

In both the dot and droplet regimes, TPST relies on 
the coherent transfer of fermionic occupation from reg- 
ister L to R. In order for this to be well-deflned, we 
choose > and Ql^Qr < so that the effective 
Dirac fermions, are conserved. In the dot regime, 
TPST can be understood by tuning to be resonant 
with a single edge mode, ^, with energy e^; so long as 
the coupling strength is weak enough to energetically re- 
solve this mode, evolution is governed by the effective 
Hamiltonian, 



eff 



^SlQIJl^u - ^9RQl/R^~k + ^-c. (6) 

and hence, state transfer proceeds via resonant fermion 
tunneling, as depicted in Fig. 4a (for details see Supple- 




FIG. 3. Fermionic spectrum of the ground state flux 
sector I The model is placed on a cylinder of circumference 
61 and width 40 unit cells, with a zigzag edge oriented such 
that the y direction is periodic [28] . The chiral edge modes are 
clearly visible near kya — n dX energies below the bulk fermion 
gap Ab = 0.46k. Numerical simulations also indicate the 
following values of the two vortex gaps: 0.1 4 At (do decagonal 
vortex) and 0.17/^ (triangular vortex). The inset illustrates 
the oppositely propagating chiral edge mode (red and green) 
at each end of the cylinder. In a droplet, the flux sector with 
completely opposite flux in each plaquette would support edge 
modes traveling in a reversed orientation. 



mentary Information) [21]. The timescale, r, required to 
achieve high fidelity state transfer depends only on the 
energy spacing between adjacent modes, Ae ^ nji^ where 
i is the system's linear dimension; to prevent the leak- 
age of quantum information into off-resonant fermionic 
modes, T>l/i<i [21]. 

In the droplet regime, we encode the fermionic occupa- 
tion into the presence/absence of a wavepacket traveling 
along the chiral edge, as illustrated in Fig. 4b [29]. Upon 
tuning both spin-registers to an energy A^, the encoding 
can be performed by choosing ^l(^) with the following 
time-dependence. 



9L{t) = 



(7) 



where f{t) characterizes the shape of the desired wave- 
packet and V is the group velocity of the chiral mode. 
Subsequent retrieval can be similarly achieved by em- 
ploying time-reversal symmetry to appropriate choose 
the shaping of QRit) (see Supplementary Information for 
details). We note that such wavepacket encoding is in 
direct analogy to the storage and retrieval of photonic 
wavepackets [30-32]. In contrast to the dot regime, the 
magnitude of the coupling strengths may be of order A^, 
which is independent of I. However, the time scale of 
TPST includes the wavepacket 's propagation time, which 



5 



depends on both the physical separation of the registers 
and the wavepacket group velocity. 



EFFECTS OF IMPERFECTIONS, DISORDER 
AND DECOHERENCE 

Having explicitly demonstrated TP ST in an exactly 
solvable model, we now consider additional imperfec- 
tions, disorder, temperature and decoherence. As the 
CSLB phase has a bulk gap and a topological invariant 
protecting its chiral edge mode, we expect the effective 
low energy fermion dynamics to be insensitive to small 
perturbations [22]. Furthermore, the chirality of the edge 
mode prevents localization and the Major ana nature of 
the edge fermions strongly suppresses the phase space for 
scattering, thereby limiting nonlinear corrections to the 
dispersion [20, 21, 24-26]. In the following, we consider 
various classes of imperfections arising from local spin 
perturbations and coupling to a finite temperature bath; 
these result in: 1) vortex excitations, 2) finite Majorana 
lifetime and 3) dynamical decoherence. 




dot regime droplet regime 



FIG. 4. Regimes of TP ST | (a) Schematic representation 
of the dot regime wherein TPST becomes analogous to tun- 
neling. In this mesoscopic dot regime, the coupling strength is 
kept weak enough to enable resolution of single edge modes, 
(b) Schematic representation of the droplet regime wherein 
TPST is achieved by mapping the quantum information from 
a spin-register onto a traveling fermionic wave-packet. The 
wave-packet is caught at the remote register, after which a 
two-qubit gate (see Supplementary Information for details) 
is performed before the information is returned (via a wave- 
packet) to the initial register. 



Finite Majorana Lifetime 



Vortex Excitations 

At low temperatures T, there will be a dilute gas of 
bulk vortices, Ny ~ n^e"^^/^, where rip represents the 
total number of bulk plaquettes. Since a vortex excita- 
tion corresponds to a tt fiux relative to the ground state, a 
circumambulating fermion acquires an additional phase 
of NyTT. Thus, the presence of vortices can have two 
relevant effects: 1) vortices within a localization length, 
^ ^ a (where a is the lattice spacing), of the edge can 
scatter a traveling fermion and 2) an odd number of vor- 
tices induces a 7r-shift of the net phase [33]. 

In addition to introducing bulk vortices, perturbations 
also generically lift the aforementioned degeneracy asso- 
ciated with dangling edge vortices. However, this will 
only affect the fidelity of TPST at the injection point, 
where one must ensure the existence of a single dangling 
edge Majorana. Away from the injection point, three 
possibilities arise: First, zero energy dangling vortices 
are completely decoupled from the fermions and hence 
will be irrelevant for TPST. Second, low energy dan- 
gling vortices will scatter only minimally, since the in- 
teraction strength between the dangling Major anas, and 
hence the hopping strength across the dangling link, will 
be extremely weak (see Supplementary Information for 
details) . Finally, much as in the bulk, the effect of high 
energy dangling vortices will be suppressed by their gap. 

As static effects, all of the aforementioned error con- 
tributions can be abrogated by the use of tomography; 
hence, it is crucial to effectively freeze out vortex fiuctu- 
ations on the time scale of TPST, and this is most easily 
accomplished at temperatures which are small compared 
to A^. 



Next, we consider the addition of generic perturbative 
local spin interactions, Hp^ to the full Hamiltonian, Ht- 
Certain classes of perturbations leave the model exactly 
solvable; more generally however, if Hp is longer ranged 
or does not respect the model's interaction symmetry, the 
gauge field acquires dynamics and the effective fermionic 
theory is no longer free. In order to understand these 
effects, we turn to a low-energy continuum theory of the 
Majorana edge (assuming that dangling vortex excita- 
tions are either decoupled or gapped out), 

He = v J ^pclcp = v j dxj{x){id)j{x), (8) 

where cj = c_p is the subset of {c^} in equation (3) 
which creates an edge excitation at momentum p and 
where we have switched to a continuum normalization of 
the Majorana field, {7(^),7(^)} = S{x — y) [33]. 

The introduction of interactions induces decay of the 
quasiparticle excitations cj. This quasiparticle lifetime 
limits the size of the droplet around which coherent ex- 
citations may be sent. The leading order symmetry- 
allowed interaction is of the form [34] 

K = >^J dx-f{x){id)-f{x){idf-f{x){idf-f{x), (9) 

where A characterizes the strength of the interaction. We 
estimate the decay rate F^"^^ of a single quasiparticle ex- 
citation using Fermi's golden rule (see Supplementary 
Information for details). In the low temperature limit 
{cp > ksT), 

F-^ ^ + + O(T^). (10) 
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To relate V^^* to the microscopic model parameters, we 
consider generic vortex-inducing local spin perturbations 
of strength which yield A ^ n{^)'^a^ in second order 
perturbation theory. Substituting into equation (10) al- 
lows us to re-express the zero temperature decay rate as 
T'^^ - ^{^Y{apY^, where A5 = vp is the energy of 
the injected TPST fermion. The surprisingly strong de- 
pendence on momenta suggests that quasiparticle decay 
can safely be neglected so long as p < 1/a. 



Dynamical Decoherence 

Finally, we consider dynamical decoherence due to 
weak coupling with a low temperature phonon bath, 
which induces additional decay T^^^ of the fermion in- 
volved in TPST. We assume that the bath couples to 
local spin operators erf and that its effect is character- 
ized by its noise spectral density, F^^'^ ^ 8(00) [35]. In 
the bulk, each such operator creates a pair of vortices 
(Fig. 2b) in addition to creating or destroying a Ma- 
jor ana quasiparticle. As the fermionic edge modes are 
exponentially localized, the contribution of this process 
to the decay rate is suppressed by e~^l^ ^ where d is the 
distance from site i to the edge. Moreover, there is an 
additional energy suppression from ^'(cJo) ~ ^-ujo/UbT 
where uoq = 2Ay is the energy cost of creating a pair of 
vortices. 

This brings us to the primary decoherence effect: edge 
noise. There are two types of spin operators acting on 
the lattice edge: 1) those which only create or destroy an 
edge fermion (type I) and 2) those which also introduce 
vortices (type H-IV), as shown in Fig. 5. Type I spin 
operators can only induce decay if they directly annihi- 
late the injected TPST edge fermion, a process costing 
energy A^. By contrast, once a vortex is created at any 
edge plaquette it can scatter the traveling TPST fermion, 
implying that the associated decoherence is enhanced by 
a factor of ^, as depicted in Fig. 5. Thus, the total TPST 
decay rate induced by edge noise is. 



-Av/ksT 



(11) 

Strikingly, the sources of decoherence in TPST are ex- 
ponentially suppressed in temperature and thus can be 
controlled [21]. 

The above analysis generalizes to other types of noise 
sources. Of particular relevance in the context of solid- 
state spin systems are nuclear spin baths, in which 
where tc is the bath's correlation time. 



In this model, the Arrhenius-type energy suppressions of 
equation (11) becomes T^^ - l/A^ if A^ > 1/tc [35]. 



EXPERIMENTAL REALIZATIONS AND 
OUTLOOK 

The search for novel topological phases represents one 
of the most exciting challenges in many-body physics; in- 
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FIG. 5. Schematic representation of the various forms 
of edge decoherence | Type I spin operators correspond to 
non-vortex inducing decoherence and can affect TPST only 
by annihilating the TPST fermion (teal star), a process sup- 
pressed by e~^*^/^^^, where As represents the detuning of 
the spin registers and hence also the energy of the injected 
quasiparticle. Type II-IV spin operators create vortices (cir- 
cular arrow) in edge plaquettes. Once a vortex is created 
at any edge site, it can scatter the traveling TPST fermion, 
leading to the decoherence being enhanced by i (the droplet 's 
linear dimension). In addition to creating vortices. Type II- 
IV spin operators also create un-gapped edge fermions (gold 
star), which we assume does not affect TPST since quasipar- 
ticle interactions have been shown to be extremely weak. 



deed, this challenge has led to a widespread effort to ex- 
perimentally identify or engineer systems exhibiting ex- 
otic topological order. One of the prototypes of such 
order is provided by the CSLB phase of the Kitaev hon- 
eycomb model; while such chiral spin liquid phases have 
yet to be experimentally implemented, several realistic 
approaches toward their realization have been envisioned. 

The realization of a honeycomb lattice, an essential 
component of implementing the Kitaev gapped B phase, 
is currently being considered in systems ranging from 
ultra-cold atoms [36-38] and polar molecules [39, 40] 
to superconducting lattices [41, 42] and dipolar-coupled 
electronic spin arrays [43, 44]. While engineering a 
macroscopic honeycomb droplet remains a daunting chal- 
lenge, recent experiments have demonstrated the ability 
to control mesoscopic ensembles containing tens of qubits 
[45-49]. Despite these remarkable advances, such meso- 
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scopic systems are insufficient in size to support the ex- 
istence of several well- separated quasiparticles, a crucial 
prerequisite to demonstrate the non-abelian braiding es- 
sential for topological quantum computing [11, 12, 22]. 
However, these smaller systems represent ideal candi- 
dates to demonstrate topologically protected state trans- 
fer and hence, the existence of a chiral fermion edge - 
another hallmark of the CSLB phase. 

Moreover, our proposed technique can also be used to 
directly characterize spin liquid states via passive spec- 
troscopy of the droplet edge. By observing the splitting- 
dependent relaxation of the spin-qubit probe, one could 
map the energy spacing between the chiral edge- modes. 
In addition, asymmetries in correlation measurements 
provide a direct indication of chirality. In this case, by 
gradually altering the physical distance separating two 
spin-qubit probes, one could characterize the timescale 
of incoherent interactions between the remote registers. 
Asymmetry in this timescale, dependent on the direction 
in which the qubits are separated provides a strong indi- 
cation of the existence of a chiral edge and would enable 
direct evaluation of the velocity associated with the edge 
dispersion. Alternatively, one could also imagine hold- 
ing the spin qubits fixed and characterizing asymmetries 
associated with L-to-R versus R-to-L TP ST. These con- 
siderations imply that solid-state magnetic spin probes 
can provide a potential tool for exploring the properties 
of natural spin liquid candidates in both organic and in- 
organic insulators [23, 50-52]. 

Finally, our technique suggests a new avenue for a 
hybrid solid-state quantum computing architecture. In 
particular, while solid-state spin systems enable the re- 
alization of long-coherence-time quantum memories and 



high-fidelity quantum registers, scaling these individual 
components up to a local area quantum network remains 
a critical challenge [14-16, 21, 43, 53]. Thus, we envi- 
sion a hybrid architecture in which conventional solid- 
state spin qubits are connected by topologically protected 
channels. In this scenario, dynamical decoupling of spins 
within an engineered CSLB droplet generates a lattice of 
mesoscopic CSLB islands around which individual qubits 
reside [54]. The edge- modes of these island regions act 
as quantum routers, ferrying quantum information be- 
tween remote spin registers. Since spin-qubits are nat- 
urally well-separated in such a hybrid architecture, the 
individual addressing and direct control of single qubits 
is greatly simplified [55]. Furthermore, such an architec- 
ture suggests a perspective in which control fields can 
reshape CSLB islands and thereby dynamically reconfig- 
ure network connectivities. 
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^ — ' I. GAUGE PROJECTION TO THE SPIN SUBSPACE 

o 

(S| Having demonstrated topologically protected state transfer by working in the extended Hilbert space in the main text, here, we 
^ consider the gauge projection back to the physical subspace. We ultimately illustrate the SWAP gate associated with TPST 
O in the language of physical spin states. Recall that the 4-Majorana representation of the spin algebra a" = 27^7^ (for spins 
a G {x^y^z} and correspondingly for Majoranas a G {1,2,3}) only holds under the constraint that D = 7^7272^0 = 1 
[SI]. This defines the two dimensional physical subspace of the four dimensional Hilbert space associated to four Majorana 
fermions. More generally, for TV spins erf, the 2^ dimensional physical space is defined by the TV constraints = 1 in the 2^^ 
dimensional extended Hilbert space. We may impose these constraints by using the physical (gauge) projector 

^ p=n^ (SI) 

^ which simply annihilates any state not satisfying the local constraint. 

The gauge projector P has several important properties. It commutes with any physical operator built out of spin operators: 
[P^af] = 0. Furthermore, P absorbs gauge transformations Di, with PDi = DiP = P. This gives us the freedom to 
understand the evolution of physical states \ip) under any physical operation O, which is a function of spin operators erf, by 
^ lifting this operation to the extended Hilbert space = (9(z7f 7q ), applying any gauge transformations which may simplify 
00 the analysis, and projecting back only at the end. That is, if we can find a state in the extended Hilbert space such that 
00 \ip) = P \ip^) then O \ip) = PO^ In particular, if \ip^) is an eigenstate of O^, the physical state \ip) is an eigenstate of 
O. This is particularly useful for identifying energy and spin eigenstates in the extended Hilbert space. One must always ensure 
^ that P ^ 0. 

o 

A. Decorated Honeycomb Lattice 

^ We now focus on the decorated honeycomb lattice model along with two additional spin registers as in Fig. 2b. of the main text. 
*^ The full Hamiltonian, Ht = Hq -\- Hl + Hr + Hint is composed of 

x,x' y-,y' z,z' 

links links links 

As ^ As 

Hl/R = -^^L 2~^^' ^^^^ ^ 9L(^L^a + 9R^R^h (S2) 

where we have chosen units in which k = 1. The extension of the Hamiltonian to the Majorana Hilbert space results in a model 
of Majorana fermions 7^^ coupled to a static Z2 gauge field (Ji^j residing on the lattice links [S2]: 

H',=\Y.^i,lh'i (S3) 

Hl,R = -^nhl-^hh% (S4) 
Hi, = -igLUL,alhl - i9RUR,blhb (S5) 
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FIG. SI: Dangling Majoranas and Vortices | A schematic microscopic edge of the decorated honeycomb with the injection point expUcitly 
labeled. Other pairs of dangling edge spins are paired into decoupled dangling vortices. It is these dangling vortices which lead to the large 
degeneracy found in the exactly solvable model. 



where Uij = i7f 7^ and a = x^y^z is the link type of (ij), or Uij = if i and j are not connected. We extend the definition of 
the gauge field Uij to the paired dangling edge spins on the boundary as in Fig. SI, and to U l,r = ^Il^r- With this choice of 
pairings, all Uij are conserved by the total Hamiltonian and thus time evolution may be understood in each Uij sector. We label 
the sectors of Uij by field configurations {Uij = ±1}. Due to the antisymmetry of Uij, there is some subtlety in correctly 
labeling sectors: in all our formulae, we take ij to be oriented according to the arrows in Fig. 2b of the main text. Thus, Uij = 1 
corresponds to a ground state gauge sector. Finally, we define Cl/r = ^ (7^/^ + ^^l/r) main text so that we may think 

of measuring the occupation of left and right register fermions. 

The Uij are gauge dependent quantities as {Di^ Uij} = 0, but the net flux around any closed loop w{C) = IlijGC ^hj E^^E^ 
invariant; thus w{C) is physical and conserved [SI]. The extended Hilbert space may be divided into conserved gauge sectors 
while the physical Hilbert space splits into conserved flux sectors after projection. As usual, we say that any plaquette P such 
that w{dP) = — 1 contains a vortex; here, we additionally extend this definition of vortices to include the dangling plaquettes 
defined by the Uij links between dangling edges, as shown in Fig. SI. These dangling vortices are completely decoupled from 
the fermions and lead to a large degeneracy of the model with open boundaries (2^^/^ where Ne is the number of dangling edge 
spins). 



B. Spin states on the Decorated Honeycomb 



Let us consider the physical ground state of the system in the absence of interaction g between the registers and the decorated 
honeycomb. The spin registers both point up, disentangled from the rest of the system, while the lattice spins sit in their collective 
ground state: It)^, \GS)q |t)^. We seek a reference ground state \ ft) in a fixed gauge sector of the extended Hilbert space such 
that \GS)q |t)^ = P 1^) up to normalization. We choose Uij = +1 (i.e. the flux configuration of the ground state sector 
as described in the main text) and we choose \ ft) to be annihilated by c^, cr and for /c > of Hq{U) + 
state is, by construction, a lowest energy eigenstate of the system, but it may not survive projection. The norm of P |^^) is 
{ft \ PP \ Q) = ^ {ft\ 1 -\- Yl - Di \ ft) where we have exploited the orthogonality of states with different gauge configurations 
Uij. The product over all gauge transformations Di measures the product of all Uij and the parity of the 7i 72 ' • • 7Ar fermionic 
state. Thus, by flipping the choice of Ul,r in |(^), we may guarantee that P \ survives projection. In general, in any fixed 
gauge sector related to our reference sector by an even (odd) number of flipped Uij's, P will annihilate states with odd (even) 
fermionic parity. 

We now construct an explicit representation of the 4 possible register states coupled to the intermediate ground state \GS)q by 
acting with cr^/^: 

lt)JG5)o|t>^ = P|J^> (S6) 

\GS)o \t)n = -iPchi \^) 
IDl \GS)o \l)n = -iPcWn \^) 
IDl \GS)o \i)n = -iPclci \^) 
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where in the last Hne we have used Ul,r = 1 acting on More generally, degenerate and/or low energy states may be 
found in either the same flux sector with (pairs of) extra edge fermions, e.g. Pc^pCp, \Q)orm degenerate flux sectors (containing 

dangling vortices), e.g. Pc^j^ \ ft) where 7f is a dangling edge Majorana. We note that since Hq{U) does not depend on the 
dangling edge Ui^s, neither do the fermionic eigenmodes Ck in these degenerate sectors nor the fermionic vacuum. 



C. SWAP Gate in the Physical Subspace 

Let us now consider time evolution U{t) in the presence of the coupling H^^^. The gauge field Uij remains conserved and the 
time evolution of the Majorana field 7^^ within each gauge sector is that of noninteracting fermions. The full Hamiltonian in our 
chosen ground state gauge sector is given by equation (5) of the main text. In general, ek, Qk,a^ and Ck depend on the gauge 
field but not on the dangling pieces of it, so the following analysis applies identically in each sector containing dangling vortices 
so long as the gauge is chosen the same way in the bulk and on UL,a^ Ur^I). Assuming that g <C A5, we may use the secular 
approximation to eliminate the c-fermion number non-conserving terms in equation (5), 

H^U) « J2 ^'^(^c^ - ^) + ^s{c{cL -\) + A5(4cH - \) (S7) 

k>0 

k 

/c k 

This Hamiltonian leaves the c-fermion vacuum \Q) invariant and evolves the modes as usual non-interacting Dirac fermions: 
it)4i2it) ' ' ' 4^(t) 1^)' where k^{t) denotes the time evolved wavefunction of the ki mode according to 
the single particle Schrodinger equation. We note that the most general Majorana evolution would mix the and c modes and 
accordingly the instantaneous c-vacuum would evolve in time. 

To enable state transfer, we now tune = for an edge mode k. In the dot regime, we further require IqlQj^ \9rQJz 6 1 ^ 
— This condition enables single-mode resolution of the edge eigenmodes and state transfer proceeds by resonant 

fermionic tunneling in an effective three mode model (dropping constants and the uninvolved modes) [S3]: 

Heff = -^dLQla'^L^k " ^dRQl/RCj, + h.C. (S8) 

Since the individual quantum registers are fully controllable, we tune ql and qr to ensure that the effective tunneling rate 
h = \ %Q~k,a\ = \ ^Q~k,h\ between the modes is equivalent. Re-expressing -iQ^ ^ = e*^^'«|(5^^| and -iQ- ^ = e'^^^'^IQ^ J, 
subsequent evolution under Hef f for a time r = results in mode evolution, 

v2tfc 

c[ -e-'*4 (S9) 

where 4> = 4>k a ~ ^k b- Using these relations to evolve the states from equation (S6), we find 

\GS)o n)n \GS% n)n (SIO) 

n)L\GS),\i)j,^ie''^\i)^\GS\n)^ 
\i}L\GS}, \i}u^-\i)L\GS%\i)^. 

up to known dynamical phases. The time evolution presented in equation (SIO) generates our desired SWAP gate in addition to 
a controlled phase gate between the register modes (up to single qubit rotations). Here \GS')q indicates a state which evolves 
from \GS)q independent of the state of the two register qubits. As depicted in Fig. S2, in combination with intra-register 
manipulations, the gate described by equation (SIO) enables universal computation between the memory qubits of the remote 
spin registers [S3]. 
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FIG. S2: Schematic representation of TPST and a remote CNOT gate | (a) The quantum registers L and R each contain two spins with 
the gold spin corresponding to the memory qubit. A single step of evolution SWAPs the information between the green spins of the left and 
right register. However, in addition to the SWAP gate, it also creates entanglement in the form of a controlled phase gate between these spins. 
To perform a remote CNOT gate between the memory qubits, we first perform an intra-register operation to SWAP the quantum information 
between the green and gold qubit of the left spin, (b) Next, the first step of TPST is performed (corresponding to a SWAP gate and a 
controlled phase gate). Afterwards, an intra-register CNOT gate between the green and gold qubit of the right register is performed, (c) The 
second step of TPST is then performed to return the information to the left register. Finally, intra-register gates are performed to yield a 
remote CNOT between the memory qubits. This enables universal computation. 



This schematic evolution holds identically for any initial state of the intermediate system \GS)q containing extra fermions or 
dangling vortices, since such states may be represented in a gauge sector where all bulk Ui^j = 1. Furthermore, in flux sectors 
in which there are an even number of bulk vortices, it is possible to choose a gauge in which Ui^j = 1 for all links near the edge. 
The evolution proceeds nearly identically in this case as well. On the other hand, in flux sectors where there are an odd number 
of bulk vortices, the energy of the edge modes is shifted by ~ k,/L implying that the spin registers are off-resonant. This can be 
corrected for through tomography and subsequent retuning [84] . 



II. SHAPING THE TRAVELING FERMIONIC WAVEPACKET 



Here, we describe the shaping of the fermionic wavepacket in the droplet regime of TPST. The edge mode energies e/^ of a 
finite-sized droplet are split at order l/£. As discussed in the main text, since single mode energy resolution becomes impossible 
in the macroscopic limit, we encode the spin register's quantum information into a fermionic wavepacket traveling along the 
chiral edge of the 2D droplet. This requires the shaping of ^l(^) and gR{t) in order to ensure the sending and receiving of the 
packet. Let us first consider the shaping of the initial wavepacket at the left register, so gR{t) = 0. As described in Sec. IC, here, 
it is sufficient for us to consider the single particle problem since the modes evolve as usual non-interacting Dirac fermions. 
By tuning A 5 to an energy in the middle of the edge dispersion and restricting \gL\ <C A5, we have (assuming a plane wave 
description of the low energy chiral edge modes) 

H^p = J2Ek\k){k\ + ^ Y.^\k){L\ + \L){k\), (Sll) 

k k 

where is the edge mode with momentum k, we have absorbed all numerical factors into gL and = vk is shifted by 
(here, we have correspondingly shifted the definition of zero energy and the indexing of k to begin at the state with energy 
As). We choose this notation for the Hamiltonian to be consistent with the literature regarding photonic wavepacket storage 
and retrieval, where an analogous problem is solved [S5, S6, S7]; thus, in this section, q, rather than being fermionic operators. 
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will represent the amplitude of the \i) mode. Initially, we consider a state {ip) whose amplitude is fully localized on the left 
spin register, \ip) = cl\L) -\- ^Ck\k), where CL(t = 0) = 1 and Ck{t = 0) = 0. After making a continuum approximation 
in both position and momentum, we formally solve the Schrodinger equation to obtain cl(^) = —^\gL{t)\'^CL(t), yielding 

cl(^) = e^^^^ where h{t) = ^ dt'\gL{t')\'^ . Substituting this result into the formal solution of Ck{t) yields 



2v JO 



Ck{t) = -i 

Thus, the shape of the outgoing wavepacket is 



i [ dt'e-'^''^*-''^^gL{t')e-''^''\ (S12) 
Jo vZ 



rt -| 

= -i / dt'5{x - v{t - t'))gL{t')e-^^''^ = -i-QL^t - x/v)e-^^'-''^^^0{t - x/v). 
Jo ^ 



where 6 is the Heaviside step function and we have assumed linear dispersion with group velocity v. Here, we note that in 
converting from a k sum to an integral, we have assumed that the amplitude on both A: < and bulk modes will be negligible 
since | ^ A5. As previously discussed, this assumption is crucial to ensure that the vacuum does not undergo time evolution. 

It is natural to think of the wave-packet in the time domain and evaluate c(x, t) dXx = 0. Thus, the solution to the problem of 
shaping any desired wavepacket, /(t), simplifies to deriving the requisite ^l(^) control function that satisfies \gL{t)^~^^^^ = 
f{t) where h{t) = ^ df\gL{t')\'^', such a solution then yields, 

SLit) = , ^^^^^ ■ (S14) 

H\\2 



The subsequent retrieval of the wave-packet at the location of the right spin register can be understood by using time-reversal; 
indeed, the control function gR{t) should be the time-reversed form of the control used to generate the time-reversed form of the 
sent wavepacket [S7]. While, for simplicity, we have considered gLj9R ^ ^ above, generalizing to complex gL/R can easily 
be achieved, for example, by employing a A-configuration spin register [S7]. One possible implementation of such a register 
is naturally provided by defect centers in diamond; in particular. Nitrogen- Vacancy color centers, which represent promising 
high-fidelity two qubit quantum registers, harbor a spin triplet (S = 1) electronic ground state [S8, S9]. 



III. THE EDGE — INJECTION AND INTERACTIONS 



There are two types of excitations on the edge: Majorana fermions and Z2 vortices [S2]. The topology of the bulk guarantees 
the existence of gapless chiral Majorana edge-modes, as described in a low energy theory by, 

He = v J dx-f{x){id)-f{x) (S31) 

where j{x) is a continuum Majorana field and v is the group velocity [SIO]. Generically all vortices are gapped; however, 
details of the lattice edge can lead to the existence of decoupled and/or low energy vortices. While the Hamiltonian (S31) does 
not capture these degrees of freedom, the presence of these additional states in the low energy Hilbert space cannot be ignored. 
Indeed, the degeneracy of the exactly solved model follows from the presence of zero-energy dangling vortices formed by pairs 
of dangling Majoranas, as depicted in Fig. SI. Away from the injection point: 1) zero energy vortices are decoupled and hence 
irrelevant to TPST, 2) low energy vortices scatter only weakly, and 3) high energy vortices are suppressed by temperature. The 
presence of a low-energy vortex degree of freedom at the injection point is critical to enable spin-edge coupling, which occurs 
at a dangling spin. Crucially, this dangling spin contains a decoupled Majorana operator 7^ (dangling Majorana), as shown in 
Fig. S 1 . Keeping track of this mode in the low-energy Hilbert space allows us to couple as follows, 

f As 

Hc = v dxj{x){id)j{x) —a^ + igcr^ J {0)^ decoupled- (S32) 

This is the continuum formulation of the microscopic coupling illustrated in the main text. To further elucidate the importance of 
the vortex degree of freedom at the injection point, we consider three possibilities. First, in the case when the dangling injection 
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Majorana is completely decoupled, the injection vortex (corresponding to the flipped Uij at the injection point) is zero energy 
and the procedure for TPST remains identical. Second, in the case when the injection Majorana is weakly interacting with a 
single nearby Majorana (respecting the interaction symmetry), the injection vortex is low-energy. In this case, the splitting of 
the spin-register will need to be retuned to account for the creation of this low-energy vortex and TPST will then naturally create 
both a fermion and an injection vortex. Crucially, tunneling of the injection vortex into the bulk will be energetically disallowed 
since /^^j^j^^t^on ^ A^^^^. Problems only arise in the third case, when the injection Majorana is interacting strongly (order 
k) with a single nearby Majorana, and the injection vortex is hence high-energy. In this case, spin-edge coupling will create 
an injection vortex which can diffuse into the bulk; thus, upon the return of the traveling fermion, the injection vortex may no 
longer be localized near the injection point, causing dephasing when the quantum information is recaptured. 



Next, we consider the role of interactions between edge modes. The fidelity of topological state transfer will be dependent 
on these interactions as they induce decay of the Majorana quasiparticles. Here, we begin by estimating the lifetime of such 
excitations in the continuum edge setting by taking into account the leading order interaction term. We consider, at T = 0, the 
situation where we tunnel a single quasiparticle excitation into the chiral edge from an associated spin register. The Hamiltonian 
is, 

H = He + H'^ (S33) 

K = ^J dx-i{x){id)-f{x){idf-f{x){idf-f{x). 

To evaluate the interaction induced decay rate of the quasiparticles, we use Fermi's golden rule and consider the relevant in- 
teraction matrix elements coupling an inconung excitation 7p|l^) with three outgoing (decayed) excitations 7pi7p2 7p3 1^)- The 
associated decay rate takes the form, 

^ = ^""Z ^I("l^-Pi'^-P^^-P3^e7pl^)l'<5(ep - e,, - e,, - e,,) (S34) 

where the delta-function imposes energy conservation and the integrals are one dimensional because the quasiparticle is confined 
to the droplet edge. To evaluate the decay rate, we begin by considering the interaction matrix element, 

nyr ^ f dki dk2 dks dk^ ^ .21'^/^ 1 1 1 \/ \ /r^^^x 

M = X ~^^^^^^^sh{^'7r)d{ki + + + A:4)(7-pi7-p2 7-p3 7/ci7/c2 7/c3 7/c4 7p) (S35) 

where we have represented the interaction Hamiltonian in momentum space; by employing Wick's theorem, we can contract 
the matrix element into a function of two point Majorana correlators. The only terms from this contraction which contribute 
are connected terms of the form {j-piJki) {j-p2lk2) il-pslks) {'yk4Jp) and such terms yield, F ~ X'^p^^/v. The leading order 
temperature correction in the limit \vp\ kT, is obtained from a Sommerfeld type expansion [Sll] and yields, 

F ^ ^ + + O(T^). (S36) 
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